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High and cold cirrus clouds mainly contain irregular ice crystals, such as, columns, 
hexagonal plates, bullet rosettes, and dendrites, and have different impacts on the climate 
system than low-level clouds, such as stratus, stratocumulus, and cumulus. The radiative 
effects of cirrus clouds on the current and future climate depend strongly on cirrus cloud 
microphysical properties including ice water content (IWC) and ice crystal sizes, which 
are mostly an unknown aspect of cirrus clouds. Because of the natural complexity of 
cirrus clouds and their high locations, it is a challenging task to get them accurately by 
both remote sensing and in situ sampling. This study presents a new method to remotely 
sense cirrus microphysical properties by using ice Raman scatter with a Raman lidar. The 
intensity of Raman scattering is fundamentally proportional to the number of molecules 
involved. Therefore, ice Raman scattering signal provides a more direct way to measure 
IWC than other remote sensing methods. Case studies show that this method has the 
potential to provide essential information of cirrus microphysical properties to study 
cloud physical processes in cirrus clouds. 
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Abstract 

To improve our understanding of cirrus cloud radiative impact on the current and future 
climate, we need improved knowledge of cirrus cloud microphysical properties. 
However, long-term studies of cirrus clouds indicated that it is a challenging task to get 
them accurately by both remote sensing and in situ sampling. This study presents a new 
method to remotely sense cirrus microphysical properties by using ice Raman scatter with 
a Raman lidar. The intensity of Raman scattering is fundamentally proportional to the 
number of molecules involved. Therefore, ice Raman scattering signal provides a more 
direct way to measure ice water content. Case studies show that this method has the 
potential to provide essential information of cirrus microphysical properties to study 
cloud physical processes in cirrus clouds. 

Introduction 

Cirrus clouds affect the surface and top-of-atmosphere energy budgets and can produce 
large local variations in atmospheric heating by heating at cloud bottom and cooling at 
cloud top. The degree and extent of the so-called greenhouse-versus-albedo effects 
involving cirrus clouds will lead to significant atmospheric differential cooling and 
heating in the vertical as well as on horizontal scales (Liou, 1986) and is dependent on 
cirrus microphysical properties and their vertical distribution (Stephens et al., 1990). 
However, it is a challenging task to measure cirrus Ice Water Content (IWC) and particle 
size by remote sensing or in situ sampling. The main limitations for in situ sampling are 
the small sampling volume and the high cost of collecting long-term data over different 
regions. IWC estimated from in situ particle size probes has large uncertainties associated 
with different ice crystal shapes and densities (Heymsfield et al. 2002). There have been 
significant advances in ground-based remote sensing of cirrus clouds using the 
Department of Energy Atmospheric Radiation Measurement (ARM) program Cloud and 
Radiation Testbed (CART) site measurements (Matrosov et al.,1994; Mace et al. 1998; 
Wang and Sassen, 2002). Nevertheless, uncertainty in the retrieved IWC by using lidar, 
radar, and radiometer measurements might be very large under some situations because 
one has to make several critical assumptions about cirrus clouds such as the size 
distribution and density of ice crystals, which vary a lot in cirrus clouds. A remote 
sensing method that possesses a signal that is directly related to IWC would be more 
promise for IWC measurements and would also be a great asset in studying these 
sensitivities in other techniques for obtaining IWC. Here, we report on a new method to 
remotely measure IWC using Raman scattering from ice in cirrus clouds. First, we briefly 
describe our Raman lidar system. Then we present the method and measurements. 
Results from a lidar-radar algorithm are used to calibrate the ice Raman scatter based 
retrieval algorithm. 


Where C is a constant and I JVC is ice water content (C IWC = c A N a ). To account for the 
possible influences of the particle morphology, the morphology factor Mis introduced. M 
is a function of the orientation, the size and the shape of the ice crystals for a given laser 
frequency. In lidar measurements, the Raman signal is integrated over a large volume 
temporally and spatially; therefore, the sensitivity of Raman scatter to ice crystal 
orientation and shape is reduced. Furthermore, Raman scattering intensity is 
fundamentally proportional to the number of molecules in the probed volume. Several 
experimental studies indicated that the Raman scatter from droplets with size parameters 
(2n r/A) larger than 60 is roughly proportional to the volume of droplets if averaged over 
a small parameter range. Therefore, we regard M as a constant for the first order 
approximation. 


Based on above discussion, the lidar equation of ice Raman scatter can be expressed as 
K 

where P ice is return power at the ice Raman scatter frequency v ict , P 0 is transmitted laser 
power, and R is the range from lidar station. C ice is a new constant including lidar system 
constant and constants of C and M in Eq. (2). or(^o,r) and a(u ice , r) represent atmospheric 
extinction coefficients including the scatter and absorption of molecules and particles at 
frequencies Vo and v !Ce and range r, respectively. 


For the Raman scatter of nitrogen, we have 


Pn,(R) = 


Cy, P 0 Nff 2 (R)cr Ni (v 0 , v N ^ , T) 
R 2 


exp(-Jo* («O 0 . r) + 0C{V ^ , r))dr). 


(4) 


where N N2 is the number density of nitrogen at range R and t/ N2 is the Raman scatter 
frequency of nitrogen. 


Combining Eqs. (3) and (4), we have 


IWC= 
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On the right side of this equation, P ice (R)/Pm(R) is known from lidar measurements; 
(C N2 N N2 (R) cr N2 (i/ 0 , VN2, T) !{C tce (Ticeiuo, Vice, T)) can be approximated as K(T) /?(R), 
where K(T) is a new constant and /?(R) is backscattering coefficient of molecules; the 
exponential term represents the difference of atmospheric attenuation difference at the 
two Raman channels and can be estimated from the Raman lidar measurements and 
atmospheric temperature profile. K might slightly depend on temperature, but we neglect 
the temperature dependency in this study. At wavelengths of v\ ct and vm, the attenuation 
difference due to clouds is negligible because of small wavelength dependency of 
attenuation for large particles comparing with wavelength. Based on these discussions, 
Eq.(5) can be expressed as: 


IWC = ^^ K p{R) & x V {^{a{v N2 ,r)-a{v ice ,r))dr), 
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( 6 ) 


Therefore, we only need to know K to determine IWC from SRL measurements, and an 
approach to determine K is presented below. 

With retrieved IWC and estimated extinction coefficient from Raman lidar 
measurements, general effective size (D ge ) can be estimated by using relationship of a 
-IWC (-2.93599* 10"*+2.54540/D ge ) developed by Fu (1996). When K is known, we 
obtain IWC and D ge profiles from SRL measurements. 

Measurements 

During the Atmospheric Infrared Sounder (AIRS) Water Vapor Experiment-Ground 
(AWEX-G) between 27 October and 16 November 2003, the SRL was sited at the ARM 
Southern Great Plains (SGP) CART site. An example of SRL signals (10-min average) 
measured on 12 November is presented in Figure 2. During this time, cirrus cloud was 
between 6 and 9 km above ground level, and optical depth estimated from elastic and 
nitrogen channels is ~1.2. A strong signal was measured in the liquid/solid Raman 
channel in the cirrus layer as indicated by strong backscatter signal and moderate 
attenuation in the A) channel. Compared with the other signals, the Raman scattering 
from water vapor and ice is weak, but, the signal-to-noise ratios of them are good for 
analysis with 10-minute average. 

As discussed above, we need to know the value of K to estimate IWC. If we know Raman 
scattering cross-section of ice and nitrogen, the optical transmittance of receiving system, 
and other system constants, we can calculate K directly. We hope to pursue this work in 
the future. Here, however, we use an alternative method to determine K by comparing 
with IWC or ice water path (IWP) derived from a published lidar-radar algorithm (Wang 
and Sassen 2002). The lidar-radar algorithm uses cirrus extinction coefficient derived 
from Raman lidar measurements and radar reflectivity factor to estimate IWC and D ge . In 
this algorithm, the accuracy of IWC or IWP mainly depends on the accuracy of extinction 
coefficient or optical depth. Between extinction coefficient and optical depth derived 
from Raman lidar measurements. Optical depth has better accuracy due to its lower 
susceptibility to multiple scattering effects (Whiteman et al., 2001). Therefore, IWP is 
better than IWC for the intercomparison to calculate the constant K. 

The AWEX-G field campaign provided a chance to use this approach to estimate K, and 
an example is given in Fig. 3. Figures 3a and 3b present time-height display of lidar 
scattering ratio from SRL and radar reflectivity factor (Ze) from the millimeter cloud 
radar at the SGP CART Site. To provide cirrus extinction coefficients for the lidar-radar 
algorithm, we simply multiply cirrus backscattering coefficients with a layer mean 
extinction-to-backscattering ratio estimated from optical depth and layer-integrated 
backscattering coefficient (Whiteman et al., 2001). IWC profiles retrieved from the lidar- 
radar algorithm are given in Fig. 3c. As mentioned above. We use IWP from the lidar- 
radar algorithm to determine K, and a comparison of IWP time series from the lidar-radar 
algorithm and from ice Raman signal is given in the Fig. 3e, which indicates that they 
agree well and follow the same pattern of variation. IWC profiles from ice Raman signal 


are presented in Fig. 3d. Comparing Figs. 3c and 3e, we see good agreement between the 
two approaches. However, there are some differences between them, which may be 
caused by the uncertainties in the lidar-radar retrieval and the lower signal-to-noise ratio 
of Raman signal at some altitudes. 

To examine the potential variation of K due to temperature and other factors, we perform 
the same calibration procedure on five different cirrus events during the AWEX-G. The 
results are listed in Table 1. These cases cover a mid-cloud temperature range of -24 to - 
60 degree with optical depth ranging from 0.05 to 2. The K calculated from each 
individual event range from 1.813 to 3.284. Comparing with the other days, the K on 2 
November 2003 has relatively low value. This is an optically thin cirrus cloud with low 
Ze, and there are potential larger errors in K because of the relatively low signal-to-noise 
ice Raman signal and relative large percentage errors in estimated optical depth and 
measured Ze. Excluding this case, K has a mean of 2.875 with a standard deviation of 
0.339, which is ~ 11.8% of the mean. Considering uncertainties in the IWP estimated 
from the lidar-radar algorithm, the variation range of K is reasonable and might mainly 
be caused by these uncertainties. Therefore, we use the mean for K. to calculate IWC. 

Applying the mean K to all cases, we can get IWC and Dge profiles. Comparisons of 
between IWP derived by the lidar-radar algorithm and that from ice Raman scattering are 
given in Fig. 4. In general, there are good agreements between these two methods. The 
intensity of Raman scattering is fundamentally proportional to the number of molecules 
involved, though the shape of the spectrum is slightly depend on temperature (Whiteman 
2003). The change of the transmitted ice Raman intensity within 40 degree is expected to 
be less than 5 %. Therefore, it is reasonable to assume that K is the same during the 
AWEX-G because of the same optical configuration. Then, the periods when the two 
IWPs do not follow each other well might indicate potential errors in the lidar-radar 
method because of possible failures of assumptions in the algorithm or errors in the 
inputs. This figure shows a potential to use Raman IWC measurements to refine other 
remote sensing methods (Matrosov et al.,1994; Mace et al. 1998; Wang and Sassen, 
2002 ). 

Conclusions 

This paper presented a new approach to remotely sense IWC and D ge using the ice Raman 
signal observed by a Raman lidar. As demonstrated in the case studies, this method 
provides the essential information of cirrus microphysical properties to study cloud 
physical processes in cirrus clouds. The main limitation of this method is weak ice 
Raman signal. But, the signal presented in this paper can be improved a lot. First, it can 
be improved by optimal design of the filter to include main ice Raman scattering 
spectrum. In the current SRL, the fdter was optimized for liquid water measurements. 
Second, it can be improved by increasing transmitted laser power. The laser power for 
measurements presented in this paper is only about 7 W. With improvements such as 
these, increases in signal to noise of between 2-4 can be expected. Therefore, Raman 
lidar might be an effective tool to measure IWC and D ge for cirrus cloud study during 
nighttime. A Raman lidar operated from an aircraft, such as the NASA/GSFC Raman 


Airborne Spectroscopic Liar (RASL at http://ramanlidar.gsfc.nasa.gov), might provide an 
effective alternative method for in situ IWC measurements with much large sampling 
volume than other approaches. Moreover, this method might provide measurements to 
refine other techniques for IWC measurements. 
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Figure Caption: 


Figure 1: Water and ice Raman scatter spectrum and the transmission of solid/liquid and 
water vapor filters used in SRL. 

Figure 2: An example of SRL signals (10-min average) measured on 12 November 2003 
at the SGP CART site. 

Figure 3: Time-height display of lidar scattering ratio (a), radar reflectivity factor (Ze, b), 
IWC from the lidar-radar algorithm(c), IWC using ice Raman signal (d), and the 
intercomparison of IWP derived from these two methods (e) observed at the SGP CART 
Site on 31 October 2003. 

Figure 4: Comparisons of IWPs derived by the lidar-radar algorithm and from ice Raman 
scatter using a mean K for the five cirrus cloud cases summarized in the Table 1. 


Table 1: The statistics of cirrus macrophysical properties and the K for the five cirrus 
events during the AWEX. 
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